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Abstract 

The  goal  of  this  paper  is  to  present  an  accurate,  comprehensive,  and  efficient,  off-road  tire  model 
for  soft  soil  applications  (traction,  handling,  ride,  and  vehicle  durability)  as  needed  to  support  cur¬ 
rent  Army  simulation  needs.  The  literature  review  revealed  that,  while  FEM  lead  to  the  most  de¬ 
tailed  tire- soil  models,  their  complexity  and  extensive  computational  effort  make  them  less  than 
ideal  for  the  applications  envisioned.  The  proposed  approach  is  a  detailed  semi-analytical  tire  model 
for  soft  soil  that  utilizes  tire  construction  details  which  parallel  commercially  available  on-road  tire 
models.  The  novelty  relies  on  increasing  the  level  of  details  for  the  tire  model,  in  improving  the  tire- 
soil  interface  model  by  enhancing  the  resolution  of  the  tire  model  at  the  contact  patch,  and  by  ac¬ 
counting  for  effects  and  phenomena  not  considered  in  existing  models.  The  model  will  be  validated 
against  experimental  data.  For  low  speed,  testing  will  be  done  on  a  single  tire  on  sandy  loam  in  the 
terramechanics  rig  at  Virginia  Tech.  The  influence  of  tire  and  vehicle  parameters  on  the  contact 
patch  forces  and  on  the  forces  transmitted  to  the  axle  will  be  investigated.  The  effect  of  soil  charac¬ 
teristics  on  the  tire  dynamics  will  be  studied.  Validation  against  data  collected  from  full  vehicle 
testing  is  included  in  the  proposed  future  work. 
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1  Introduction 

The  goal  of  this  paper  is  to  present  an  ac¬ 
curate,  comprehensive,  and  efficient,  off-road 
tire  model  for  soft  soil  applications  (traction, 
handling,  ride,  and  vehicle  durability)  as 
needed  to  support  current  Army  mobility 
goals.  The  literature  review  revealed  that, 
while  finite  element  method  (FEM)  leads  to 
the  most  detailed  tire-soil  models,  their  com¬ 
plexity  and  extensive  computational  effort 
make  them  less  than  ideal  for  the  applications 
envisioned.  The  proposed  approach  is  a  de¬ 
tailed  semi-analytical  tire  model  for  soft  soil 
that  utilizes  tire  construction  details  which 
parallel  commercially  available  on-road  tire 
models.  The  novelty  relies  on  increasing  the 
level  of  details  for  the  tire  model,  in  improv¬ 


ing  the  tire-soil  interface  model  by  enhancing 
the  resolution  of  the  tire  model  at  the  contact 
patch  and  by  accounting  for  effects  and  phe¬ 
nomena  not  considered  in  existing  models. 
The  model  will  be  validated  against  experi¬ 
mental  data.  For  low  speed,  testing  is  con¬ 
ducted  on  a  single  tire  on  sandy  loam  in  the 
terramechanics  rig  at  Virginia  Tech.  The  in¬ 
fluence  of  tire  and  vehicle  parameters  on  the 
contact  patch  forces  and  on  the  forces  trans¬ 
mitted  to  the  axle  will  be  investigated.  The 
effect  of  soil  characteristics  on  the  tire  dynam¬ 
ics  will  be  studied. 

2  Literature  Review 

Most  tire  models  are  developed  with  a 
particular  application  in  mind.  Depending  on 
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the  application,  a  compromise  needs  to  be 
made  between  computational  efficiency,  accu¬ 
racy,  resolution,  and  ease  of  parametrization. 
Fig.  1  illustrates  the  classification  of  tire  mod¬ 
els  based  on  the  approach  used,  in  increasing 
order  of  complexity,  and  their  applications.  At 
the  bottom  of  the  diagram  are  the  least  com¬ 
plex  models  such  as  the  empirical  Pacejka- 
type  models;  at  the  very  top  are  the  very  com¬ 
plex  finite  element  models. 
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Fig.  1  :  Tire  model  types  and  their  intended  applications  [1] 

A  brief  overview  of  the  tire  models  in 
each  of  these  classes  is  presented  next. 

2.1  Empirical  Tire  Modelling 

Empirical  tire  models  rely  on  modelling 
the  tire  performance  by  curve-fitting  a  very 
large  amount  of  experimental  data.  Most  of 
the  time  the  mathematical  formulations  thus 
obtained  do  not  carry  a  physical  significance, 
and  are  specific  for  each  tire  simulated.  The 
most  important  areas  of  application  for  empir¬ 
ical  models  are  mobility  and  real-time  simula¬ 
tions.  Perhaps  the  most  important  empirical 
tire  model  is  the  Waterways  Experimental 
Station  (WES)  mobility  model  [2].  This  is  a 
ugo/no  go ”  model  that  uses  an  instrument 
called  a  cone  penetrometer  to  determine  a 
cone  index,  which  is  then  compared  to  a  mo¬ 
bility  index  of  a  vehicle.  The  vehicle  mobility 
index  is  computed  based  on  vehicle  character¬ 
istics  such  as  vehicle  weight,  contact  area, 
size  of  grouser,  engine  power,  and  type  of 
transmission.  If  the  vehicle  mobility  index 
exceeds  the  cone  index  then  motion  is  possi¬ 
ble;  otherwise  the  vehicle  is  stuck. 

Wismur-Luth  [3-4]  also  developed  an 
empirical  model  using  bias-ply  tires.  When 


radial  tires  are  used,  the  Wismur-Luth  model 
tends  to  under-predict  the  traction  of  the  vehi¬ 
cle.  Most  models  that  originated  from  the 
Wismur-Luth  model  used  variations  of  the 
Wismur-Luth  equations  to  account  for  chang¬ 
es  in  experimental  conditions,  including  the 
use  of  radial-ply  tires.  Brixius  [5]  developed  a 
traction  prediction  model  for  bias -ply  tires 
based  on  similar  methods,  and  suggested 
modifications  of  the  relations  for  radial  tires. 
However,  Upadhyaya  and  Wulfsohn  stressed 
in  [6]  that  the  empirical  traction  models  sel¬ 
dom  provided  any  insight  into  the  underlying 
mechanics  and  should  be  used  with  caution 
when  evaluating  new  tire  situations. 

2.2  Semi-analytical  Tire  Modelling 

Semi-analytical  tire  models  are  probably 
the  most  versatile  models  because  they  pro¬ 
vide  good  fidelity,  while  also  having  physical 
significance.  Their  large  spectrum  of  com¬ 
plexity  makes  them  suitable  for  a  large  variety 
of  applications.  The  simplest  semi-analytical 
models  are  the  single  point  contact  models, 
mostly  used  for  ride  applications.  There  are 
also  complex  semi-analytical  tire  models  used 
for  a  variety  of  vehicle  dynamics  simulations. 
Most  semi-analytical  tire  models  simulations 
run  close  to  real  time. 

The  FTire  [7]  is  arguably  the  most  widely 
used  commercial  semi-analytical  tire  model 
for  on-road  applications.  This  model  discre¬ 
tizes  the  tire  in  N  discrete  masses  (in  a  single 
plane)  to  model  the  tire.  By  using  torsional 
springs  and  shape  functions  the  model  pre¬ 
dicts  the  tire  behaviour  in  the  lateral  direction. 
A  simple  contact  model  based  on  Bekker’s 
equations  is  available  for  off-road  simulations. 
Other  tire  models  based  on  a  similar  semi- 
analytical  approach  are  the  CD-Tire  40  [8] 
and  RMOD-K  [9]. 

Chan  [10]  has  previously  developed  flex¬ 
ible  ring  and  rigid  tire  models  that  are  used  for 
on-road  and  off-road  vehicle  dynamics  and 
vehicle  mobility  applications.  Lee  et  al.  [11] 
developed  a  semi-analytical  model  for  snow 
where  the  snow  is  modelled  as  a  depth- 
dependant  modified  Drucker-Prager  material. 
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This  model,  as  well  as  Chan’s  model,  use  a 
combined  slip  formulation.  El-Gawwad  et  al. 
[12-13]  developed  a  multi-spoke  tire  model 
using  Bekker’s  pres  sure- sinkage  formula  and 
included  the  effect  of  straight  lugs  on  vehicle- 
terrain  interaction  for  an  agricultural  tire.  This 
tire  model  accounts  for  parameters  which  af¬ 
fect  vehicle  maneuvering,  such  as  slip  angle, 
soil  deformation  modulus,  lug  dimensions, 
and  lug  spacing.  Harnisch  developed  a  model 
for  both  rigid  and  elastic  tires,  where  the  elas¬ 
tic  tire  model  “ uses  a  larger  substitute  circle 
to  describe  the  deformed  contact  patch  be¬ 
tween  the  tyre  and  soiF  [14].  The  substitute 
circle  diameter  is  obtained  from  the  equilibri¬ 
um  between  the  vertical  reaction  force  of  the 
soil  and  the  normal  force  applied  to  the  tire. 

2.3  Finite  Element  Tire  Modelling 

The  use  of  the  finite  element  methods  to 
model  tires  flourished  in  the  last  two  decades, 
mostly  because  of  the  increase  in  computing 
power  available.  Fervers  [11]  created  a  model 
that  uses  a  multi-layered  ring  to  model  the 
different  components  in  the  tire.  Given  the 
complexity  and  time  required  to  create  and 
run  a  three  dimensional  model  he  transfers  3D 
elements  into  a  2D  half-space  to  save  on  com¬ 
puting  time.  The  results  from  this  study  are 
encouraging  in  the  sense  that  it  demonstrates 
that  a  detailed  model  that  includes  the  differ¬ 
ent  layers  of  the  tire  adds  to  the  resolution  of 
the  model;  proving  that  a  detailed  model  can 
capture  with  a  good  resolution  the  tire  dynam¬ 
ics  of  off-road  simulations. 

Nakashima  and  Oido  [12]  discussed  a 
method  to  combine  Finite  Element  Analysis 
(FEA)  and  Discrete  Element  Analysis  (DEA) 
to  model  tire-soil  behavior.  The  soil  is  mod¬ 
eled  using  DE  while  the  tire  is  modeled  using 
FE.  The  main  purpose  of  this  study  is  to  better 
describe  the  microscopic  deformation  of  the 
soil.  Lee  et  al.  [30]  also  employs  an  FE  model 
for  fresh,  low-density  snow.  In  this  study  the 
snow  is  simulated  as  a  pressure-sensitive 
Drucker-Prager  material  in  order  to  establish 
the  relationship  between  strengths  in  com¬ 
pression/tension. 


One  of  the  most  representative  FE  tire 
models  is  Shoop’s  model  [15].  For  simula¬ 
tions,  a  Drucker-Prager  snow  model  is  used. 
Composite  material  properties  are  used  to 
model  the  different  parts  of  the  tire.  The  mod¬ 
el  predicts  the  forces  at  the  wheel  axle,  but  it 
is  still  computationally  expensive. 

2.4  Literature  Review  for  Tire  Testing 

Throughout  the  years  of  studies  in  vehicle 
performance,  particularly  on  deformable  soils, 
a  variety  of  instrumented  vehicles,  trailers, 
and  indoor  single-tire  testing  platforms  have 
been  constructed  to  aid  in  understanding  the 
tire-soil  interaction.  Using  such  test  facilities 
and  equipment  it  was  possible  to  gather  rele¬ 
vant  test  results.  It  is  important  to  be  familiar 
with  these  past  studies  in  order  to  design  fu¬ 
ture  experiments.  Since  the  Advanced  Vehicle 
Dynamics  Lab  at  Virginia  Tech  has  a  single 
wheel  testing  unit,  most  of  the  testing  meth¬ 
ods  reviewed  here  will  revolve  around  that  of 
a  single  wheel  test. 

The  overall  focus  of  building  a  controlled 
tire  testing  device  is  to  investigate  tire  per¬ 
formance,  such  as  traction,  under  known  nor¬ 
mal  loads  and  controlled  slip  conditions.  The 
information  collected  supports  tire  models 
validation  and  helps  understand  the  tire  per¬ 
formance  under  given  conditions.  It  may  also 
help  in  controlling  soil  compaction,  improv¬ 
ing  the  protection  of  the  environment,  increas¬ 
ing  soldier  safety,  and  optimizing  the  tire  de¬ 
sign  for  performance  and  energy  efficiency. 

Among  all  the  different  types  of  tire  test¬ 
ing  platforms,  indoor  soil  bins  appear  to  be  an 
appropriate  experimental  method  to  follow, 
especially  for  terramechanics,  due  to  the  ca¬ 
pability  of  controlling  the  soil  conditions. 
This  proved  necessary  in  early  field  tractor 
tire  testers  which  would  suffer  from  soil  dis¬ 
continuities  [16].  A  single  wheel  test  device 
consists  of  a  carriage  that  can  move  at  its  own 
speed  simulating  a  vehicle,  via  a  motor-chain 
system,  and  houses  the  test  wheel  and  any  sort 
of  instrumentation  of  interest.  The  carriage 
can  also  house  the  motor  that  drives  the  wheel 
independently;  this  allows  the  possibility  of 
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testing  any  desired  longitudinal  slip.  The  car¬ 
riage,  which  acts  like  a  quarter-car  model 
without  a  suspension  system,  can  run  a  test 
tire  over  the  conditioned  soil  in  the  soil  bin. 

The  indoor  soil  bins  allow  the  investigator 
to  mix  and  compact  the  soil  as  desired,  run  a 
test,  then  repeat  the  reconditioning  process. 
However,  constantly  reconditioning  the  soil 
after  every  test  run  and  having  limited  space 
can  make  this  method  tedious  and  constrict¬ 
ing.  The  mix  and  compaction  procedure 
would  begin  by  tilling  the  soil,  leveling  it,  and 
then  applying  a  load  to  compact  the  soil  to  a 
desired  uniform  density.  Rather  than  manual¬ 
ly  reconditioning  the  soil  with  rakes  and  roll¬ 
ers,  some  studies  have  developed  automated 
devices  that  process  the  soil  as  needed  before 
each  test  run  as  done  in  the  study  at  Kyoto 
University  in  Japan  [17].  In  this  agricultural 
study,  a  device  was  fitted  in  the  single  wheel 
test  rig  to  till,  level,  and  compact  the  water- 
filtration  type  sand.  This  “Mixing-and- 
Compaction”  had  a  tiller  in  the  front  and  an 
aluminum  plate  in  the  rear  which  would  flat¬ 
ten  the  sand  surface  [17].  To  operate  the  de¬ 
vice,  a  motor-driven  roller  was  attached  in  the 
rear,  and  it  was  also  used  to  compact  the  soil. 

The  type  of  tire  test  would  also  dictate 
how  the  soil  should  be  reconditioned,  as  dis¬ 
cussed  in  the  study  from  the  Department  of 
Biological  and  Agricultural  Engineering  in 
University  Putra  Malaysia,  UPM  [18].  In  this 
study,  two  different  soil  reconditioning  meth¬ 
ods  were  used  with  either  a  driven-wheel  or  a 
towed-wheel.  Both  compaction  procedures 
began  by  loosening  up  the  sand  to  at  least  150 
mm  depth  with  a  hand  hoe,  then  leveling  it  by 
hand  [19].  If  in  driven- wheel  mode,  steel- 
tubed  rollers  filled  with  concrete  were  hand 
rolled  across  the  soil  bin  [18].  If  in  towed- 
wheel  mode,  a  heavy  gasoline-powered  soil 
compactor  was  run  across  the  soil  surface 
equipped  with  a  vibrating  plate. 

Aside  from  reconditioning  the  soil  after 
every  test  run,  some  tests  have  been  per¬ 
formed  to  disregard  that  step  in  order  to  ex¬ 
plore  multi-pass  effects  [20].  The  goal  of  this 
sort  of  test  is  to  simulate  the  rear  wheel  of  a 


vehicle  rolling  along  the  same  path,  or  rut,  left 
by  the  front  wheel. 

Performing  tests  such  as  a  multi-pass  tests, 
involves  several  factors;  one  key  factor  that 
the  investigator  can  control  is  the  longitudinal 
speed  of  the  wheel  carriage.  The  longitudinal 
slip  typically  ranges  from  0-100%,  but  can 
also  go  to  -100%  [2].  When  performing  con¬ 
stant-slip  tests,  the  carriage  and  test  wheel  go 
at  their  respective  velocities  at  a  constant  rate, 
as  defined  by  the  investigator  at  the  start  of  a 
test,  for  the  effective  rolling  radius.  Although 
the  drawbar  pull,  Fd,  defined  as  the  difference 
between  the  tractive  effort  and  the  resultant 
resistance  force,  can  be  controlled  as  well  in 
certain  single  wheel  devices,  it  was  concluded 
from  a  Waterways  Experiment  Station  (WES) 
study  that  controlled-slip  tests  are  better  than 
controlled-pull  tests  for  defining  pull-slip  and 
torque-slip  relations  [16,2].  It  is  also  interest¬ 
ing  to  note  that  in  this  WES  study,  wheel  per¬ 
formance,  as  described  by  controlled-slip 
tests,  was  found  to  be  independent  of  the  start¬ 
ing  conditions  once  the  selected  slip  is  at¬ 
tained  [16].  This  conclusion  has  been  rein¬ 
forced  by  the  recent  study  conducted  by 
Woodward  at  NASA  GRC  [21], 

From  all  the  studies  reviewed,  including 
some  that  are  not  discussed  here  due  to  space 
restrictions,  it  is  well  understood  that  using  a 
systematic  methodology,  with  known  influ¬ 
encing  factors,  it  is  possible  to  correlate  single 
wheel  testing  results  to  those  obtained  from 
full- vehicle  testing  [16]. 

3  Model  Overview 

The  proposed  soft  soil  tire  model  is  in¬ 
tended  for  mobility,  traction,  ride,  handling, 
and  durability  applications,  for  which  a  com¬ 
plex  semi-analytical  tire  modelling  approach 
is  the  best  fit,  when  also  taking  into  account 
its  computational  efficiency.  The  schematic  of 
the  proposed  model  can  be  observed  in  Fig.  2. 
It  is  a  discretized  lumped  mass  approach  that 
uses  springs  and  dampers  in  multiple  configu¬ 
rations  to  represent  the  different  sections  of 
the  tire.  The  model  is  structured  in  three  paral¬ 
lel  planes,  two  of  them  representing  the  side- 
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walls,  and  one  of  them  represents  the  belt  and 
the  tread  of  the  tire.  By  differentiating  be¬ 
tween  the  sidewall  and  the  belt  a  more  realis¬ 
tic  application  of  the  local  forces  can  be  im¬ 
plemented.  Each  plane  consists  of  Nm  number 
of  masses;  Nm  is  suggested  between  80  and 
100,  but  it  is  user  defined.  This  approach  has 
the  advantage  of  being  modular  and  it  allows 
the  user  flexibility  in  determining  the  resolu¬ 
tion  of  the  model. 

Each  lumped  mass  has  three  degrees  of 
motion:  the  translational  motion  in  all  direc¬ 
tions.  Moreover,  there  is  relative  motion  be¬ 
tween  masses  in  the  same  plane  and  between 
masses  in  adjacent  planes,  as  well.  Relative 
motion  is  also  allowed  between  lumped  mass¬ 
es  and  the  rigid  wheel  in  the  circumferential 
direction,  which  is  a  novelty  introduced  by 
this  model.  The  wheel  has  all  six  degrees  of 
freedom.  The  total  degrees  of  freedom  of  the 
model  are, 


DOF  =  9  Nm  +  6  (i) 

The  shape  of  the  contact  patch  is  as¬ 
sumed  to  be  a  curved  trapezoid  to  accommo¬ 
date  transient  situations  such  as  sudden  steer¬ 
ing  maneuvers;  during  quasi-steady-state  situ¬ 
ations  the  shape  will  switch  to  a  curved  rec¬ 
tangle,  which  according  to  the  literature  [19- 
21]  is  a  good  assumption  for  such  conditions. 


Fig.  2:  3-Dimensional  view  of  the  tire  model 


3.1  Sidewall  Element 

The  sidewall  element  is  more  complex 
than  the  belt  and  tread  element  because  it  is 
directly  connected  to  the  rigid  wheel.  A  dia¬ 
gram  of  the  model  can  be  observed  in  Fig.  3. 
The  mass  of  a  single  sidewall  element  is  cal¬ 
culated  using  Eq.(2). 

0.5  rnsidewall  (^) 


Fig.  3:  Sidewall  diagram.  Left-  the  in-plane  connections; 
right  -  the  out-of-plane  connections 


3.2  Belt  and  Tread  Element 

The  belt  and  tread  element  is  similar  to 
the  sidewall  element,  except  it  lacks  the  con¬ 
nection  to  the  rigid  wheel.  The  diagram  for  the 
belt  and  tread  element  can  be  found  in  Fig.  4. 


Fig.  4  :  Belt  and  tread  element  diagram.  Left  -  the  in-plane 
connections;  right  -  the  out-of-plane  connections. 


Proceedings  of  thel7,k  International  Conference  of  the  ISTVS  -  Sept.  18-22,  2011.  Blacksburg,  VA,  USA 


5 


UNCLASSIFIED 


It  is  also  important  to  note  that  the  belt 
and  tread  plane  has  a  larger  radius,  thus,  ac¬ 
counting  for  the  curved  shape  of  the  tire. 

The  mass  of  a  single  belt  and  tread  ele¬ 
ment  is  calculated  using  the  equation, 


m 


™t,bt  = 


tread 


T  mbelt 


(3) 


Nr, 


p(0)  =  k0z(9)nz(9)m 


(4) 


p(0)  =  k0z(9)n 


14 

1  -sd 


sin 


arccos 


z0  -  z(0)V 

Rm  ) 


(5) 


4.2  Shear  Stress 


4  Tire  -  soft  soil  interaction 

4.1  Pressure-sinkage  relationship 

One  of  the  most  important  aspects  of  the 
interaction  of  the  tire  with  soft  soil  is  the 
computation  of  the  tire  sinkage.  In  this  study 
we  implement  an  approach  developed  by 
Grahn  [23],  which  is  a  method  derived  from 
Bekker’s  formulation  [2],  but  with  additional 
features;  it  includes  the  effects  of  longitudinal 
slip,  vehicle  velocity,  and  vertical  penetration 
velocity.  The  results  by  Grahn  shown  in  Fig.  5 
demonstrate  the  large  influence  that  penetra¬ 
tion  velocity  has  on  the  pressure-sinkage  rela¬ 
tionship,  thus  motivating  the  need  to  account 
for  it  in  our  model. 


The  shear  stress  in  the  contact  patch  is 
computed  using  Janosi  and  Hanamoto’s  ap¬ 
proach  [2].  In  this  approach  the  longitudinal 
shear  displacement  is  calculated  at  each  mass 
by  integrating  the  tangential  slip  velocity  over 
time  yielding  the  following  equation, 

jx(9)  =  Ri(dme  ~  9)  (6) 

-  (1  -sd)(sin  9e  -sin0)] 

On  the  other  hand,  the  lateral  shear  dis¬ 
placement  is  calculated  by  integrating  the  lat¬ 
eral  velocity  over  time,  thus,  the  equation  be¬ 
low, 

jy{9 )  =  /?i(0)(l  —  sd)(9e  —  9)  tanac  (7) 


sinkage /mm 

Fig.  5  :  Pressure-sinkage  relationships  for  different  vertical 
penetration  velocities  [2] 


Using  this  methodology,  the  sinkage  is 
computed  for  each  mass  in  the  discretized  tire 
model.  The  pressure  in  the  radial  direction  on 
each  mass  can  be  computed  using  Eq.  (4)-(5) 
Note  that  if  m  —  0  then  these  equations  yield 
Bekker’s  quasi-static  formulation.  All  the  var¬ 
iables  are  defined  in  the  nomenclature  table. 


Once  the  shear  displacements  are  calcu¬ 
lated,  the  shear  stress  can  be  evaluated  using 
Bekker’s  approach  [24]  at  each  mass  ? 

/  -M±\  (8) 

T(0)  =  Tmax  (l  -  e  k  j 

Nonetheless,  in  order  to  account  for  a 
combined  slip  scenario,  the  shear  strength 
envelope  of  the  soil  needs  to  be  defined.  This 
is  done  by  implementing  the  Mohr-Coulomb 
failure  criteria,  where  the  inequality  in  Eq.  (9) 
bounds  the  shear  stress. 

TX  \2  /  Tycp  \2  <  j. 

T-max'  sT max ' 

Umax  =  (c  +  crn(0)  tan(0))  (10) 
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Substituting  Eq.  (10)  into  Eq.  (9)  and  9a 
for  9  in  both  Eqs.  (6)  and  (7)  yields  the  fol¬ 
lowing: 


-R;[(0e-0a)-(l-sd)(sin  (9e-sin  Ga)] 


1  —  e 


+  1  —  e 


-( Ri(l-Sd)(Se-ea )  tan  ac) 


(11) 


=  1 


where  9a  is  the  angle  of  transition  be¬ 
tween  adhesion  and  sliding.  Thus,  by  solving 
numerically  for  9a  in  Eq.  (11)  the  shear  stress 
can  be  obtained  using  Eqs.  (11),  (12)  and  (13). 


The  sinkage  from  the  maximum  stress 
angle  to  the  exit  angle  is  calculated  using  Eq. 
(16), 


z(B)  =  Rl(B)[  cos  6t 


(16) 


/  6  —  6r  \ 

~  KeT^eJ {de  ~  6m) , 


Kdm  ~  Or- 


COS  0, 


The  maximum  sinkage  is  evaluated  by 
Eq.  (17). 


Zo  =  RltfmXc OS  9m  ~  COS  Be) 


(17) 


/  -R1[(0e-0a)-(l-sd)(sin0e-sin0a)]\ 

TmaxU-e  J  oe>e>ea  (12) 

T-max  @a  ->  @  —  ~^r 

/  -(Ri(l-sd)(0e-eg)  tan  ac)\ 

W  1  -  e  ky  Ge>G>Ga 

V  /  (13) 

T-max  >  G  >  ~Gr 

4.3  Normal  Stress 

For  rigid  driven  wheels  the  location  of 
maximum  normal  stress  is  located  below  the 
wheel  axle,  which  is  the  lowest  part  of  the 
tire.  Experimental  results  show  that  the  loca¬ 
tion  of  maximum  pressure  for  a  pneumatic  tire 
is  somewhere  between  the  entry  and  exit  an¬ 
gles  and  its  location  is  dependent  on  slip  ratio 
[24].  Thus,  in  order  to  account  for  this  phe¬ 
nomena,  Wong’s  approach  [24]  is  to  deter¬ 
mine  the  sinkage  of  each  mass  by  first  deter¬ 
mining  the  angle  of  maximum  normal  stress 
9m.  The  experimental  constants  Co  and  cm  are 
both  a  function  of  normal  load  and  inflation 
pressure. 


For  straight  line  driving  at  zero  camber,  a 
parabolic  distribution  of  the  pressure  in  the 
contact  patch  in  the  longitudinal  direction  is 
an  accurate  representation.  The  shape  of  the 
contact  patch  in  this  case  is  close  to  rectangu¬ 
lar.  Most  existing  tire  models  assume  a  uni¬ 
form  distribution  of  the  pressure  in  the  lateral 
direction  of  the  contact  patch.  In  this  study,  a 
parabolic  distribution  will  be  considered  in  the 
lateral  direction  during  steady- state  maneu¬ 
vers.  However,  for  cornering  maneuvers  and 
driving  with  a  non-zero  camber  angle,  a  trape¬ 
zoid  is  a  more  accurate  representation  of  the 
contact  patch;  in  this  case,  the  pressure  distri¬ 
bution  in  the  lateral  direction  is  no  longer 
symmetric,  as  it  is  assumed  for  straight  line 
driving  at  zero  camber. 

The  first  step  in  the  simulation  is  to  load 
the  tire.  Thus,  by  substituting  Eqs.  (14)-(17) 
into  Eq.  (5)  and  using  a  quasi-steady  state 
equilibrium  for  the  tread  and  belt  masses  in 
the  radial  direction,  the  following  equation  is 
obtained: 


Gm  =  [co  (Fz>Vt)  +  Cm(Fz,pt)\sd\]8e  (14) 

Based  on  this  angle  of  maximum  normal 
stress  the  sinkage  is  defined  in  two  sections. 
The  sinkage  from  the  entry  angle  to  the  max¬ 
imum  stress  angle  is  determined  from  the  fol¬ 
lowing  formulation, 

z(0)  =  /?((0)(cos  6  —  cos  0e)  (15) 


mi%i  ^6 \2Ri'i  Rl,i, sidewall, left 

R-l.i.sidewall.right )  —  ~^i  sin  0j  +  an  i  +  (18) 

U:,;sin0( 

Where  the  radial  direction  is,  in  fact,  the 
normal  direction  in  which  the  pressure  is  ap¬ 
plied  on  each  mass  element. 

vn(8)  =  P(0)  (19) 
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By  replacing  Eq.  (19)  into  (18)  an  equa¬ 
tion  with  both  Rj  and  (),,  is  obtained.  Thus, 
additional  information  is  needed,  in  order  to 
approximate  the  initial  loaded  radius  of  the 
tire.  The  formulation  used  to  approximate  the 
initial  loaded  radius  of  each  mass  is  based  on 
the  approach  developed  by  Sandu  and  Chan  in 
[16].  The  main  difference  is  that  in  the  current 
study  the  sidewalls  and  the  tread  and  belt  ele¬ 
ments  will  have  different  effective  rolling  ra¬ 
dii,  since  they  are  located  in  different  planes. 
Based  on  this  technique  it  is  possible  to  solve 
for  the  entry  angle  and  consequently  for  the 
maximum  stress  angle  and  the  corresponding 
sinkage  of  each  mass  element. 

4.4  Bulldozing  Effect 

The  bulldozing  effect  is  created  when  a 
volume  of  the  soil  in  the  contact  patch  is  dis¬ 
placed  to  the  sidewall  of  the  tire.  As  such,  a 
lateral  force  is  created  on  the  sidewall  of  the 
tire.  A  graphical  representation  of  this  phe¬ 
nomenon  can  be  observed  in  Error!  Refer¬ 
ence  source  not  found.. 

Few  tire  models  account  for  the  bulldoz¬ 
ing  component  of  the  lateral  force.  Those  that 
incorporate  it,  applied  it  at  the  single  mass 
that  represented  the  contact  patch.  In  the  cur¬ 
rent  study,  due  to  the  discretization  of  the  tire, 
the  bulldozing  effect  can  be  applied  directly 
on  the  sidewall  elements,  thus  increasing  the 
realism  of  the  model. 


Fig.  6:  Graphical  representation 


Thus,  a  bulldozing  force  will  be  applied 
at  each  mass,  which  is  also  a  more  realistic 
representation.  The  formulation  used  to  de¬ 
termine  the  bulldozing  effect  is  based  on  the 
principle  of  passive  ground  resistance  devel¬ 
oped  by  Terzaghi,  which  is  presented  in  [22]. 
This  formulation  has  different  dimensionless 
soil  resistance  coefficients  that  according  to 
Wong  [2]  are  dependent  on  the  soil  angle  of 
friction  (p.  On  the  other  hand,  q  is  the  sur¬ 
charge  load  from  the  accumulated  soil,  which 
is  calculated  by  assuming  that  the  soil  dis¬ 
placed  in  the  lateral  direction  on  each  side  of 
the  wheel  is  the  same  as  the  volume  of  soil 
accumulated  for  the  surcharge  load. 

Fybd.i  =  (YsZiNy  +  cZiNc  +  qiZiNq)  cos  Sf  <2°) 

_  Ys^soil  (29 

^ 1  mz(0) 


4.5  Multi-pass  Effect 

The  multi-pass  is  a  really  important  effect 
for  soft  soil  interaction.  The  impact  of  a  tire 
on  soil  is  dependent  on  a  variety  of  factors. 
One  such  factor  is  the  type  of  pass.  A  towed 
wheel  will  not  alter  the  properties  of  the  ter¬ 
rain  in  the  same  manner  as  a  driven  wheel, 
which  induces  larger  changes  [24-25]. 


of  the  bulldozing  effect  and  shear  stresses 
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The  number  of  passes  also  has  an  influ¬ 
ence  on  the  terrain  properties.  Moreover,  the 
slip  ratio  also  affects  the  soil  properties  for 
multi-pass. 

The  approach  used  in  this  study  to  quan¬ 
tify  the  multi-pass  effect  will  rely  on  the  work 
of  Holm  [26]  and  Senatore  [27],  and  its  main 
idea  is  to  estimate  new  parameters  for  the  soil 
based  on  the  pass  type  and  slip. 

4.6  Forces  and  Moments 

4.6.1  Longitudinal  Force/Drawbar  Pull 

The  drawbar  pull  is  one  of  the  most  im¬ 
portant  performance  metrics  in  off-road  driv¬ 
ing,  and  it  gives  the  measure  of  the  capability 
of  the  vehicle  to  perform  work,  such  as  pull¬ 
ing  a  trailer  or  a  plow.  The  drawbar  pull  is 
computed  in  the  same  way  for  all  the  masses 
in  the  three  planes  using  Eq.  (23).  The  total 
drawbar  pull  will  be  computed  using  Eq.  (24). 


3  Nm 

Fy  ~  ^  Fy:Sidewaii  i  +  Fy jjgit  i  (27) 

i= 1 

4.6.3  Driving  Torque 

The  driving  torque  is  calculated  using 
Eq.(28),  where  W,-  is  the  wheel  load  at  each 
mass. 

3  Nm 

T  =  ^  Fl,i[rx,i  ~  Wi  cos  °i\  (28) 

;= i 

4.6.4  Overturning  Moment 

The  overturning  moment  at  the  axle  is  a 
great  indicator  of  performance  in  handling 
situations  and  it  is  calculated  using  Eq.  (29). 

3  Nm 

^  —  Fy  (Ri  i  COS  (29) 

i= 1 


Fx,i  =  DPi  =  Txi  cos  Bi  -  an,t  sin  et  (23) 


3  Nm 


;= i 


4.6.2  Lateral  Force 

The  calculation  of  the  lateral  force  is  also 
very  important,  especially  for  steering  maneu¬ 
vers.  Note  that  the  lateral  force  is  made  up  of 
two  components,  which  are  produced  by  the 
(i)  lateral  shear  stress  and  by  (ii)  the  bulldoz¬ 
ing  effect.  For  the  sidewall  elements  the  lat¬ 
eral  force  is  calculated  using  Eq.  (26),  while 
for  the  belt  and  tread  the  lateral  force  is  com¬ 
puted  using  Eq.  (25).  The  total  lateral  force  is 
given  by  Eq.  (27). 

Fy,belt,i  —  TyCpyi  (25) 

Fy, sidewall, i  ~  Fybd,i  T  ~^ycp,i 


4.6.5  Self-Aligning  Torque 

The  calculation  of  the  self-aligning 
torque  at  the  axle  follows  the  same  approach 
of  the  overturning  moment,  and  is  done  using 
Eq.  (30). 

3  Nm 

'  Fn^Fy  i  sin  9 i 
;= i 


5  Model  Validation 

In  order  to  insure  the  accuracy  of  the 
proposed  tire  model  and  provide  the  semi- 
analytical  model  with  input  parameters,  it  is 
imperative  that  well-controlled  experimental 
testing  is  performed. 

The  proposed  design  of  experiments  for 
the  validation  testing  is  thoroughly  explained 
in  this  paper.  The  testing  will  be  primarily 
performed  on  a  single-wheel  testing  platform, 
the  Terramechanics  rig  in  the  Advanced  Vehi- 
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cle  Dynamics  Lab  (AVDL)  at  Virginia  Tech. 
More  detailed  information  about  the  design 
and  instrumentation  on  the  rig  are  found  in 
[28].  The  Terramechanics  rig  was  designed  to 
simulate  various  road  surfaces,  to  apply  nor¬ 
mal  load,  and  to  measure  the  forces  and  mo¬ 
ments  caused  by  the  tire-soil  interaction. 
Testing  can  also  be  performed  at  different  toe 
and  camber  angles,  as  well  as  at  different  soil 
compaction,  moisture,  and  depth  levels  to 
replicate  real-world  conditions. 

A  statistically  relevant  design  of  experi¬ 
ments  was  drafted  to  obtain  optimal  data.  The 
Terramechanics  rig  carriage,  shown  in  Fig.  7, 
allows  a  tire  to  be  tested  while  the  data  acqui¬ 
sition  program  on  the  central  computer  allows 
the  user  to  vary  certain  parameters  during  the 
run.  A  key  parameter  that  will  be  varied  in  a 
wide  range  of  values  is  the  longitudinal  slip  of 
the  tire,  which  is  controlled  by  the  two  sepa¬ 
rate  motors,  of  the  rig  carriage  and  of  the  ac¬ 
tual  wheel,  respectively. 


Fig.  7  :  Terramechanics  rig  carriage  equipped  with  a  pneu¬ 
matic  tire 


Before  testing  on  the  sandy  loam  selected 
for  this  study,  a  series  of  tests  on  a  rigid  sur¬ 
face  have  been  performed  to  insure  the  correct 
calibration  and  the  accuracy  of  the  sensors 
used  on  the  rig.  The  hardware  on  the  Ter¬ 
ramechanics  rig  is  inspected  and,  if  needed, 
calibrated  after  every  run,  in  order  to  avoid 
any  drift  in  the  data  collection. 

Figure  Fig.  8  illustrates  the  drawbar  pull 
results  obtained  recently  while  testing  a 
Goodyear  P255/55R18  104  H  generic  winter 
tire  at  a  5.18  kN  normal  load  on  the  rigid  plat¬ 
form  in  the  rig. 


_ 

^  

Raw  Data 

Fitted  Curve 

- £ 

Slip  (%) 

Fig.  8:  Normalized  drawbar  pull  plotted  against  slip 

By  inspecting  Fig.  8,  one  can  notice  that 
the  normalized  drawbar  pull-versus-slip  curve 
follows  the  expected  behavior  of  a  tire  on  a 
rigid  surface.  While  the  tests  on  sandy  loam 
are  still  to  be  conducted,  this  preliminary 
work  is  illustrative  of  the  type  of  experimental 
results  to  be  obtained,  as  well  as  of  the  large 
range  of  slip  ratios  that  can  be  employed. 

5.1  Test  Parameter  Selection 

The  design  of  experiment  is  mainly  fo¬ 
cused  on  varying  the  longitudinal  slip  and  the 
normal  load.  Because  we  are  interested  in 
observing  how  these  parameters  affect  the  tire 
performance,  different  levels  per  parameter 
will  be  tested.  As  recommended  by  the  indus¬ 
try  quad  members  of  the  project,  the  tire  be¬ 
havior  at  longitudinal  slips  larger  than  65% 
will  be  investigated  since  it  is  at  such  high  slip 
ratios  that  the  tire  performance  on  soft  soil  is 
less  understood.  The  tires  will  be  tested  at  two 
normal  loads  (100%  and  150%  of  the  nominal 
load)  and  at  two  slip  angles. 

5.2  Tires  Selection 

The  experimental  study  is  scheduled  to 
be  performed  on  two  different  types  of  tires. 
One  of  the  tires  is  the  Goodyear  Wrangler 
SilentAmor  P265/70R18  (81.3  cm  diameter) 
tire.  This  is  a  large,  all-terrain,  light-truck  tire 
that  is  highly  rated  for  its  performance  on  dry, 
wet-,  and  off-road  applications  which  rank  9, 
9,  and  8  out  of  10,  respectively,  on  the  Good- 
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year  website  [29].  The  second  is  a  Goodyear 
Rawhide  III  which  is  a  smaller  (63.5  cm  di¬ 
ameter),  all-terrain  vehicle  (ATV)  tire.  This 
tire  has  a  rugged,  deep-bite  tread  pattern  made 
for  soft  terrain  applications,  super  rigid  side- 
wall,  and  is  the  common  OEM  tire  choice  for 
many  Gators,  ARGO,  and  Honda  ATVs. 

5.3  Soil  Parameters 

The  process  of  reconditioning  the  soil  af¬ 
ter  every  test  run  will  certainly  be  strenuous 
for  the  experiment,  as  it  will  mostly  be  exe¬ 
cuted  by  hand,  so  the  minimal  amount  of  vari¬ 
ance  in  soil  parameters  is  best.  The  soil  of 
choice  in  this  study  is  sandy  loam,  which  is  a 
good  compromise  between  non-cohesive  sand 
and  very  cohesive  clay.  An  amount  of  7.26 
metric  tons  of  sandy  loam  donated  to  the  pro¬ 
ject  by  Egypt  Farms  provides  sufficient  soil 
for  testing  the  tires  of  interest  without  the  re¬ 
sults  being  affected  by  the  boundary  condi¬ 
tions  (rigid  surface  under  the  soil  and  the  rig 
side  walls). 

Two  different  soil  moisture  content  levels 
between  0%  and  30%  will  be  used. 

Another  soil  parameter  of  high  interest  is 
the  soil  compaction.  A  Rimik  cone  penetrom¬ 
eter  CP  40  II  is  employed  to  test  the  compac¬ 
tion  resistance  of  the  soil.  Using  generic  sandy 
loam  data  found  in  [5],  a  simplified  form  of 
Bekker’s  equation  was  used  to  estimate  a 
normal  range  of  compaction  resistance  for 
such  a  soil.  Within  this  range,  the  compaction 
resistance  between  333  N  and  445  N  has  been 
chosen  for  testing. 

5.4  Other  Parameters  of  Interest 

To  help  validate  tire  forces  during  steer¬ 
ing  maneuvers,  tests  must  be  performed  for 
non-zero  camber  and  toe  angles.  The  rig  was 
designed  to  accommodate  a  camber  angle 
between  0  and  8  degrees  and  a  toe  angle  from 
-45  to  45  degrees.  The  levels  of  camber  se¬ 
lected  for  this  study  are  0  and  4  degrees.  For 
the  toe,  the  levels  are  -15  and  +15  degrees. 


5.5  Initial  Design  of  Experiment 

All  the  parameters  previously  mentioned 
have  been  compiled  into  the  following  design 
of  experiments  testing  matrix  in  :  Initial  de¬ 
sign  of  experiment  test  matrix.  The  design  of 
experiments  will  be  processed  using  statistical 
software  to  insure  statistically  relevant  and 
repeatable  results.  The  software  will  provide 
the  recommended  number  of  runs  per  test 
[30].  In  deciding  the  final  number  of  runs 
conclusions  of  past  studies  conducted  at  WES 
will  also  be  considered  [2,16] 


Table  1:  Initial  design  of  experiment  test  matrix 


Parameter 

No.  of 

Levels 

Range 

Slip  Percentage 

8 

20  -  100% 

Normal  load 

2 

100%  and  150% 

nominal  load 

Inflation  Pressure 

2 

75%  and  100%  max 

pressure 

Moisture  Content 

2 

0  and  30% 

Compaction  Re¬ 
sistance 

2 

333-445  N 

Toe  Angle 

2 

0  deg  and  30  deg 

Camber  Angle 

2 

0  deg  and  4  deg 

6  Conclusions 

This  paper  presents  a  proposed  semi- 
analytical  soft  soil  tire  model  and  its  experi¬ 
mental  validation  methodology.  The  model  is 
intended  for  traction,  ride,  handling,  and  du¬ 
rability  applications.  The  work  on  the  project 
is  ongoing  and  simulation  results  as  well  as 
validation  tests  will  be  presented  in  future 
publications. 
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Nomenclature 

b  Tire  width 

Initial  radius  approximation  dimen- 
p  sionless  variable 

c  Soil  cohesion 

Experimental  constant  used  to  de- 
C°  terrnine0m 

Cj  Multipass  experimental  parameter 

csd  Multipass  experimental  parameter 

Experimental  constant  used  to  de- 

Q 

m  terrnine0m 

8  Tire  deformation 

Longitudinal  shear  deformation 
*  modulus 

ky  Lateral  shear  deformation  modulus 

kxi  Multipass  experimental  parameter 

ksd  Multipass  experimental  parameter 


c1  Sidewall  radial  damping  (in-plane) 

Sidewall  circumferential  damping  (in- 
C2  plane) 

Lateral  inter-element  damping  (out- 
°3  of-plane) 

Wheel-sidewall  radial  damping  (in- 
64  plane) 

Rh  Loaded  radius 

Ru  Un-deformed  radius 

sd  Longitudinal  slip 

14  Vehicle  longitudinal  velocity 

z0  Maximum  sinkage 

n  Pressure-sinkage  index 

m  Penetration  velocity  exponent 

m bead  Bead  mass 

mibt  Belt  plus  tread  element  mass 

m thread  Tread  mass 

mbelt  Belt  mass 

^sidewall  Total  Sidewall  mass 

mis  Sidewall  element  mass 

mwheel  Wheel  mass 

Nm  Number  of  masses 

ern  Normal  stress 

tx  Longitudinal  shear  stress 

Ty  Lateral  shear  stress 

jx  Longitudinal  shear  displacement 


k0  Static  modulus  of  soil  deformation 

Sidewall  radial  spring  stiffness  (in- 

1  plane) 

Wheel-sidewall  radial  spring  stiffness 

2  (in-plane) 

Inter-element  radial  spring  stiffness 

3  (in-plane) 

^  Tangential  inter-element  spring  stiff- 

4  ness  (in-plane) 

Lateral  inter-element  spring  stiffness 

5  (out-of-plane) 

Radial  inter-element  spring  stiffness 

6  (out-of-plane) 


jy  Lateral  shear  displacement 
p  Radial  pressure 

pt  Tire  internal  pressure 

0  Soil  angle  of  internal  friction 

y  Soil  density 

ys  Unit  weight  of  soil 

YSi  Multipass  experimental  parameter 

ysd  Multipass  experimental  parameter 
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Ny  Soil  specific  weight  coefficient 

Nc  Soil  cohesion  coefficient 

Nq  Soil  surcharge  load  coefficient 

Surcharge  load  from  accumulated 
^  bulldozed  soil 

z  Sinkage 

Initial  radius  approximation  dimen- 
’  sionless  variable 

z  Vertical  penetration  velocity 

z0  Maximum  sinkage 

Sf  Angle  of  friction  between  soil  and  tire 

Central  angle  describing  mass  posi¬ 
tion 

9a  Angle  of  transition 

9b  Trailing  Edge  Angle 

9e  Entry  angle 

9f  Leading  edge  angle 

9m  Angle  of  maximum  normal  stress 

9r  Exit  angle 

ac  Slip  angle 

y  Lateral  distance  from  tire  center 

f  Toe  Angle 
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